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ABSTRACT 

We report here the discovery of an new near-equal mass Trans-Neptunian 
Binaries (TNBs) L5c02 and the and the putative detection of a second TNB 
(L4kl2)among the year two and three detections of the Canada-France-Echpic 
Plane Survey (CFEPS). These new binaries (internal designation L4kl2 and 
L5c02) have moderate separations of 0.4" and 0.6" respectively. The follow- 
up observation confirmed the binarity of L5c02, but L4kl2 are still lack of more 
followup observations. L4kl2 has a heliocentric orbital inclination of ~ 35°, 
marking this system as having the highest heliocentric orbital inclination among 
known near-equal mass binaries. Both systems are members of the classical main 
Kuiper belt population. 

Based on the sample of objects searched we determine that the fraction of 
near-equal mass wide binaries with separations > 0.4" is 1.5% to 20% in the 
cold main classical Kuiper belt and, if our detection the binarity L4kl2 holds, 
3% to 43% in the hot main classical objects are binary. In this manuscript we 
describe our detection process, the sample of objects surveyed, our confirmation 
observations. 

Subject headings: Kupier belt. Trans- Neptunian object. Binary 



1. Introduction 



Jus t as the detection of t he first trans-Neptunian object (TNO, neglecting Pluto), 1992 
QBi, by I Jewitt fc Luu I ( 1l992l ) let to an avalanche of new TNO discoveries, the serendipitous 
discovery of the first trans-Neptunian binary (TNB, neglecting the Pluto-Charon system). 



1998 WW31, by IVeillet et al. I (|2002[ ) spurred t he discovery of do zens more such systems. 
There are currently more than 50 TNBs known (INoU et al. Il2008al ). The growing catalogue 
of known TNOs and TNBs discovered in characterized surveys permits the use of these 
objects as probes of structure formation and evolution in the outer solar system. 



The large sky separation, 1.2 and the similar brightness, A 



mag 



0.4, of the two com- 
ponents of the 1998 WW31 provided the first clue to the enigmatic nature of TNBs. That 
this was the first kind of TNB discovered, however, is partially due to the observational re- 
quirement of being able to split the two sources in-order to determine that one is observing a 
binary. Since the TNOs are rather faint, most observations are of rather low SNR (~ 10) and 
to detect a TNB in such observations generally requires that the two components have an 
angular separation of order the seeing of the observations. The next TNB discovered, 2001 
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QT297 (lOsip et al. Il2008l ) also had large separation that was at the limit of the resol ution 
of the observations. The discovery separation of 2001 QW322 (~ 4". |Petit et al. II2OO8I ) was 
much larger than the seeing of the observations, perhaps indicating that large separation 
binaries are more than just an artifact of the observing circumstances. At this time about 
~ 10% of the known TNBs have separation in excess of 1". The pathway of wide separation 
TNBs are, as yet, not well understood. In general, the study of binary systems provides 
very direct information on the physical nature of the constituent members. Binary systems 
provide the only reliable means to derive the masses of TNOs. Furthermore, if their sizes can 
be computed by using albedos obtained from optical and infrared measurements, their aver- 
age bulk densities can be estimated to to roughly 10% accuracy. While these measurements 
provide data on the physical characteristics of the binary components, knowledge of the 
physical properties for the TNO population provide inputs to understanding the formations 
of the TNB themselves. 

In the cold classical disk of TNOs inclinations < 5°, the occurrence rate of TNBs with 



(Stephens et al. 


2006; 


Noll et al. 


2008a 


b) 



be the result of the dynamical histories of different TNO populations, indicative of different 
formation histories or survival rates for TNBs in the two populations. To help understand 
the difference in TNB properties and thus unravel the formation conditions of the outer solar 
system more attention should must be given to the census of different types of TNBs (i.e., 
wide binaries vs. close binaries) in different regions of the Trans- Neptunian belt. 

Several large-scale ground-based surveys of TNOs have been performe d or in pro- 
cess. These include the Deep Ecliptic Survey (DES) at CTIO and Magellan (IMillis et al. 



2002 



Canada-Frarice E cliptic Plane Survey (CFEPS, (jjones et al. 



2006 



Elliot et al. I l2005l ) , the Keck Telescope survey jTruiillo fc Brown I boosh. and the 



Kavelaars et al. 2009 



Petit et al. 



2010h) . The first two surveys have yielded three TNBs out to a total of 362 



TNOs (INoU et al. Il2008al ). In this work we will report the putatitve detection of two new 
TNBs based on the CFEPS observations. The observations used in our TNB search and our 
search technique are described in Section 2.1 and 2.2. The results of our search and follow-up 
observations are given in Section 2.3 and 2.4. Section 3 provides a discussion and summary. 



2. Observation and Analysis 

The Canada-France Ecliptic Plane Survey (CFEPS) used the observations of the Very 
Wide sub-component of the Canada- France-Hawaii Telescope Legacy Survey (CFHT-LS). 
The CFEPS Kuiper belt study is a primary scientific goal of the CFHLS-VW (Legacy Survey- 
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Very Wide) which employs the MegaCam camera with a 1" x 1° field of view. The purpose 
is to provide a Kuiper Belt survey with quantitatively-known detection biases and long- 
term tracking to achieve precise orbits of KBOs. CFEPS has surveyed fields whose ecliptic 
latitudes range between +/ — 2° from the ecliptic plane. Approximately 500 square degrees 
have been covered with a detection limit of about 24.3 magnitudes in g' and subsequently 
re-imaged, for tracking, in r' and i'. At completion of the CFHT-LS component of the 
survey approximately 200 TNOs had been discovered a nd tracked to thre e oppositions, here 



we report on the search for TNBs among these TNOs (jPetit et al. Il2010l ). 



The CFEPS images of newly discovered TNOs can be used to determine if any of the 
TNOs are in fact TNBs. The most direct method is to look for separate components of 
the binary systems, as has been done by previous authors. A visual search of the discovery 
images for TNBs is only sensitive to large (> 1") separation TNBs, because of the limitation 
in angular resolution of the detector (~ 0.187" per pixel) and seeing condition (> 0.7"). No 
TNBs were found during the initial visual inspection of the CFEPS discovery images. 



2.1. Point Source Modeling 

To improve our sensitivity to smaller separation TNBs (< 1") we examine and com- 
pare the point-spread functions of individual KBO images to check if they are truly point 
sources or possible extended sources similar to that employed in previous TNB searches 
jKern and EUiot Ibood : Istephens et"al~ll2006h . The IRAfQ DAOPHOT package was used 



to build a point-source model (PSF) for each of the three CFHT/MegaCam images in which 
a CFEPS TNO was discovered. Several (10-15) isolated bright(SNR > 50) point-sources 
were hand selected as the model inputs. The DAOPHOT package computes the average 
shape of the point sources to obtain the PSF model (Figure [1] ). The model input stars 
were selected to be within ~ 1' from the target TNO to avoid the distortion across the the 
CFHT/MegaCam field. 

The utility of the PSF modeling approach to binary detection can be determined via 
construction of artificial binaries. The IRAF/ADDSTAR task was used to create artificial 
near equal magnitude TNBs by adding a series of closely separated doublet sources whose 
total fiuxes are similar to those of the TNOs being examined. These simulation reveal that 
TNBs with separations larger than about 80% of the FWHM of the image have goodness- 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Asso- 
ciation of Universities for Research in Astronomy (AURA) under cooperative agreement with the National 
Science Foundation. 
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of-fit values more than 2-o" larger than the for point sources in the same image. 
Artificial binaries with separations smaller than 80% of the FWHM for the image could not 
be distinguished from point sources using this method. 

For each image of a given TNO the goodness-of-fit parameter was determined. 
Sources whose value of differed significantly from that of similar brightness point sources 
in the image were flagged as possible TNBs. Visual examination of the flagged candidates 
revealed that in all but one case the observed large residuals were caused by proximity to 
background sources and not a close companion. In one case, L5c02, the deviant value 
appeared to be associated with a close companion (see Figure [2]). The large number of false 
positives and the low rate of successful detection is likely due to the combination of low 
SNR of the CFEPS TNO detection images coupled with the low-frequency of large (> 0.7") 
separation TNBs. 



2.2. Source elongation 

Examination of the image shape along a single spatial-dimension using a simple func- 
tional form, provides an improved approach for binary detection. This approach reduces the 
number of degrees of freedom (as compared to the PSF fltting approach) while increasing 
the available signal by collapsing the source proflle along one spatial dimension. A one- 
dimensional Gaussian provides a reasonable match to astronomical images and has only 
three degrees of freedom: total flux, centroid, and width. 

Using the ellipse fltting routine in IRAF/ IMEXAM the ellipticity and position an- 
gle of each TNO can be determine. The flux along the minor axis of the TNO image is 
then collapsed and the best-flt Gaussian shape parameters along the major axis are deter- 
mined. The width of the Gaussian that best matches the collapsed TNO light proflle is 
then compared to the distribution of widths for point sources in the image, also collapsed 
along the same direction as the TNO image. The deviation of the width computed for the 
TNO from the expected value based on the reference point sources is used to flag possible 
TNBs. Those sources with width more then 2-0" larger than the mean for point sources in 
the same frame were examined visually. Again, all TNOs flag for visual inspection where 
found to have background contamination except for two sources: L5c02 (previously detected 
as binary using PSF modeling) and perhaps L4kl2. Tests with artiflcially created near equal 
magnitude TNBs revealed that this one-dimensional fltting approach is sensitive enough to 
detect binaries with separations larger than 1/3 the FWHM. 
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2.3. Detected TNBs 

Figure [2] presents a comparison of the image subtraction and elongation fitting of three 
CFEPS TNOs: L4j05, L4kl2 and L5c02. Two of the sources in Figure [2] (L4j05, L4kl2) 
have residuals that are normal for point sources in their images and values of consistent 
with a point source. The TNO L4j05 is found to have an elongation width consistent with 
point sources detected in the same frame and is presented as an example of a none-detection. 
L4kl2, however, exhibits significant deviation from the corresponding one-dimensional Gaus- 
sian fitting of point sources in the image, revealing the extended nature of this source. The 
TNO L5c02 exhibits obvious residuals compared to the PSF model and also exhibits one- 
dimensional Gaussian shape which is much broader than the points sources in the same 
frame. 

Using these two procedures (PSF modeling fitting and profile elongation) we identified 
two binary candidates (L5c02 and L4kl2) out of 63 TNOs examined. Those 63 TNOs are the 
brightest TNOs from CFEPS samples (r < 23.0). Table [T] gives the orbital parameters and 



photo metric data of these two objects. Based on the classification system of iGladman et al. 



( I2OO8I ). both L5c02 and L4kl2 are members of the classical main Kuiper belt. That is, neither 
are scattering or have a resonant orbit wit h Neptune, and they h ave semimajor axes between 
the 3:2 and 2:1 mean-motion resonances (IGladman et al. II2008I ). 



2.4. Follow-up Observations 

In order to confirm the binary nature of L4kl2 and L5c02, we have obtained high- 
resolution follow-up images using the Subaru 8-meter and WIYN 3.5-meter observatories. 

Observations of L5c02 obtained with the WIYN telescope on 3 May 2008 are presented 
in Figure El An elongated structure with a length of about 0.67" (21,000 km at 44.2 AU) is 
clearly revealed in this image. The two components from this image are not fully resolved, 
however, simultaneous PSF modeling using two points sources reveals that this object is well 
represented by two sources of comparable magnitudes. 

Observations of L4kl2, obtained with Subaru on April 23, 2009, are shown in Figure Si 
An elongated structure with a major axis of about 0.8" (25,000 km at 43.6 AU) is quite 
prominent in the images. Although the 0.4" to 0.8" separation change is plausible, it is 
possible this single Subaru image is of a background star behind the TNO. If the elongated 
structure is not due to background star behind the TNO, this image is well modeled by 
doublet source with both sources having comparable magnitudes. 
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Table [2] provides a summary of the information on the CFEPS discovery, recovery and 
follow-up observations. 



3. Discussion and Summary 

The images of the new TNOs discovered by CFEPS have been examined with a view to 
detect the presence of binaries. Both the methods of PSF subtraction and a newly developed 
scheme of 2D Gaussian profile fitting have been used. We only examine the objects brighter 
than g ~ 23, because both two methods are SNR sensitive. Of the 63 new TNOs (31 of 
63 are classical belt objects, 20 are the resonant objects) examined, two, L4kl2 and L5c02 
have been identified to be wide binaries and L5c02 have been confirmed. Both of them are 
classical TNOs. 

To estimate the binary fraction in different dynamical classes, especially in cold and hot 
components of classical Kuiper belt, it is important to realize that an z < 5 degree orbit 
could be drawn from either a low-inclination or higher-inclination component when both 



are modeled as gaussians. We used the CFEPS survey simulator (iKavelaars et al. II2009I ) 
to determine that for ^ 23 about 7% of the z < 5° detections are actually from the hot 
component, and thus it is problematic to label a particular TNO as 'cold' just because it has 
z < 5 degrees. Conversely, if L4kl2 does turn out to be binary, with heliocentric i ~ 35° the 
probability it is drawn from the narrow inclination 'cold' population is essentially negligible. 
In our 31 classical belt objects sample, 23 of 31 are i < 5°. Therefore we expect that ~ 21 
are cold classical belt objects and 10 are hot components. For the 95% confidence interval, 
based on the Poisson uncertainty, we have constrained the binary rate to between 1.5% and 
20% for the 'cold belt' and between 3 and 43% for the hot belt. Also, the 95% confidence 
interval for detections is an upper limit of 3 sources, we know the binary fraction among 
the 'resonant' sources is below 3/20, or 15%. 

This fra ction is substantia lly lower than that determined from high-resolution HST 
Noll et al. I 2008"b[) but is c onsistent with the fraction of wider separation sys- 



imagmg [see 



tems reported in lKern and Elliot I (120061 ). We note that L4kl2 is has the highest heliocentric 



orbital inclination among all currently known binary sources. 

The present numerical algorithm will be a promising approach in the search of TNBs in 
large surveys like Pan-STARRS which is expected to multiply the number of detected TNOs 
many times. A systematic search as described in the present work will shed new light on the 
orbital distribution of TNBs and their origin. 



Additional high quality tracking observations of L4kl2 from GEMINI (A.H. Parker, 
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private communications) have not resolved (seeing FWHM ~ 0.5") this source as binary. 
Given the lack of a foUowup observations of our Subaru confirmation data we caution that 
L4kl2 may actually be a tight binary (separation < 0.5") discovered in the CFEPS images 
during its apocentre passage. We are continuing to monitor this source in an attempt to 
fully resolve the nature of this putative TNB. 

This work was supported in part by NSC Grant: NSC 97-21 12-M-008-011-MY3 and 
Ministry of Education under the Aim for Top University Program NCU. 
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Table 1. Heliocentric orbital data and system magnitudes 



Object 


a(AU) 


e 


i (degree) 


Dynamical Class 


Magnitude 


Mutal Period(day)* 


L4kl2 


40.77 


0.117 


35.23 


Classical 


23.0 (r') 


320 


L5c02 


45.74 


0.035 


1.791 


Classical 


22.8 (r') 


314 



Note. — *Assume that semi-major axis equal to separation of two components with 
Igjcrr? bulk density. 
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( ) r^Gl^TTrCl "I" 

Uoci VcLL 


lUli oLlllillicLi y 




Telescope 


UT Date 


A mag 


Separation(") 


PA(degree)* 


L4kl2 










CFHT 


2005 May 14.4 




0.4 


270 


Subaru 


2009 Apr 23.6 


0.71 


0.79 


270 


L5c02 










CFHT 


2005 Feb 10.3 


0.28 


0.61 


28 


CFHT 


2005 Apr 11.3 


0.26 


0.62 


35 


CFHT 


2006 Feb 06.3 


0.13 


0.75 


95 


WIYN 


2008 May 03.2 


0.32 


0.67 


207 



Note. — PA is angle east of north. 




Fig. 1. — An example of the PSF of single stars in a CFEPS field which is constructed by 
averaging the PSFs of several single stars around a TNO: (a) a cross-section of the PSF; (b) 
a two-dimensional projection. 
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Original CFEPS Image Subtraction Image PSF Profile 



(a) L4jd5 July 08, 2004 

■ ■■ ■ ^ 

. , .;■ . ® ■ ■■ 

■■. ' ■ ' . ■ 

: ■ "i " ." " ' "■ 


(d) L4j05 


^^^^ 

(b) L4kl2 May 14,- 2005 

J . ■ ■ . ' ■ ^' 

:. ^ -■ ■ ~-- . 

.- r" -■. : ; 


(e) L4ki2 . ^ ■■ .: 

■ ■ ■ ■ ■ • ■ 

J ■ ■ . - ■ =■ ■ .■ ° 
^ ■ . 1 ■ ■ 

■ . r . '■ ■ ! " - L .. 


(c) L5c02 N<^. 04, 2005 

-r.;. ■■■■■■■ -i^V^^^^ 


(f) L5c02 ^ ... . . . 

■ ■. 
• ■ ^' 1- .• - r ■ 








(i) L5c02 




Fig. 2. — A comparison of (a-c) the original images of three TNOs L4j 05 (non-detection), 
L4kl2 and L5c02, which are detected by CFEPS, (d-f) their images after subtraction from 
the reference PSF, and (g-i) the fitting of their PSF profiles(sohd hnes) along the elongated 
axis to the reference PSFs(dash lines). The normalize pixel values of TNOs are shown in 
circle dots. 
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,.;:?r^- V..^ . ....... / .V 



iib ■ ^. Il?y- 3., 200.^ . .•■Apr. ll,20q5 . ; 

Fig. 3. — The confirmation image of L5c02 taken by WIYN (3.5 m at Kitt Peak) on May 
3, 2008. Also sliown are tlie CFEPS discovery image (Apr. 1, 2005) plotted in the same 
angular size for comparison. 




Fig. 4. — The single image of L4kl2 taken by Subaru on April 23, 2009. Also shown are the 
CFEPS discovery image (May 14, 2005) plotted in the same angular size for comparison. 



